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Abstract 

State sum constructions, such as Kuperberg’s algorithm, give partition functions of 
physical systems, like lattice gauge theories, in various dimensions by associating local 
tensors or weights, to different parts of a closed triangulated manifold. Here we extend 
this construction by including matter fields to build partition functions in both two and 
three space-time dimensions. The matter fields introduces new weights to the vertices 
and they correspond to Potts spin configurations described by an A-module with an 
inner product. Performing this construction on a triangulated manifold with a boundary 
we obtain the transfer matrices which are decomposed into a product of local operators 
acting on vertices, links and plaquettes. The vertex and plaquette operators are similar 
to the ones appearing in the quantum double models (QDM) of Kitaev. The link operator 
couples the gauge and the matter fields, and it reduces to the usual interaction terms in 
known models such as Z 2 gauge theory with matter fields. The transfer matrices lead to 
Hamiltonians that are frustration-free and are exactly solvable. According to the choice 
of the initial input, that of the gauge group and a matter module, we obtain interesting 
models which have a new kind of ground state degeneracy that depends on the number 
of equivalence classes in the matter module under gauge action. Some of the models have 
confined flux excitations in the bulk which become deconfined at the surface. These edge 
modes are protected by an energy gap provided by the link operator. These properties 
also appear in “confined Walker-Wang” models which are 3D models having interesting 
surface states. Apart from the gauge excitations there are also excitations in the matter 
sector which are immobile and can be thought of as defects like in the Ising model. We 
only consider bosonic matter fields in this paper. 
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1 Introduction 


Lattice systems with gauge and matter fields is a recurring theme in high energy and condensed 
matter physics. They are a useful way of regularizing field theories in particle physics where 
the standout example is lattice QCD, a widely studied field [Tj. The gauge groups involved in 
these programs are compact, continuous non-Abelian groups like 577(3) and 51/(2) and Abelian 
groups like U( 1). Whereas in the context of condensed matter physics such lattice systems are 
used to model magnetism in solids and crystals and study various phases of strongly correlated 
electrons most of which are typically low energy phenomena. One of the prominent examples 
is the Hubbard model which is an interacting microscopic model used to describe a variety of 
strongly correlated phenomena [2]. In these examples the gauge groups considered are finite 
discrete groups such as the Abelian Z n . The classic examples using such Abelian groups as 
degrees of freedom are the Ising chain, Heisenberg model and the Potts model [3]. An important 
connection between lattice spin systems and gauge theories was given in the famous article by 
Kogut bq®- They study in detail the connection between a d-dimensional quantum spin system 
and the d + 1-classical system via the transfer matrix method. This gives the correspondence 
between the two-dimensional quantum Ising model and three dimensional Z 2 lattice gauge 
theory. 

The partition functions of such lattice systems can be found out in a systematic way by 
looking at them as state sum models. These state sum models are used in several areas of 
mathematics and physics, including statistical mechanics [6], random matrices \7 \ , knot theory 
[8j, lattice gauge theory [9] and quantum gravity |T0]. Such state sum models are defined based 
on a combinatorial decomposition of a manifold such as a lattice or a triangulation, which 
can be interpreted to be the space-time in a physical picture. In these models the degrees of 
freedom live on the vertices and/or edges of the lattice. In a more general setting one could 
associate local states on even faces and volumes. We restrict our attention to those living only 
on the vertices and edges. A state in such a system is the tensor product of the configuration 
in each local edge and vertex. Local weights are associated to the vertices edges and faces 
of the triangulated lattice. These weights can be thought of as tensors with indices which 
can be raised and lowered. These tensors are then allowed to contract and thus one obtains 
a number corresponding to the partition function when this construction is carried out for a 
closed triangulated manifold. 

These constructions have been carried out in the past for finding 3-manifold invariants which 
were also interpreted as the partition functions of certain physical models. They were especially 
useful in constructing the partition functions of Topological Quantum Field Theories (TQFTs). 
An important example is the 3D Dikgraaf-Witten invariant HU which furnished all topological 
lattice gauge theories. Hamiltonian realizations of these invariants in two dimensions is given 
in G2- Another important example is the Turaev-Viro type of TQFTs in 3D [13]. These are 
realized by the Levin-Wen models or the string-net Hamiltonians ra in two dimensions. These 
are exactly solvable models made up of commuting projectors describing degrees of freedom 
belonging to a unitary fusion category located on the edges of the lattice. They realize anyons 
as low energy excitations and provide a general class of long-ranged entangled topologically 
ordered phases in two dimensions. The Turaev-Viro invariant are known to be equivalent to 
the Chain-Mail link invariant m It was shown in [16] how one can obtain the string-net 
models from the Chain-Mail link invariant which is also a knot invariant. 

In the spirit of such computations we showed in m that the partition function of Quantum 
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Double Hamiltonians (QDM) of Kitaev [T31IT9] can be obtained by a deformation of Kuperberg’s 
3-manifold invariants 123. The QDM Hamiltonians realize lattice gauge theories based on an 
involutory Hopf algebra A. Lattice gauge theories based on a group G is a special case of this 
as A = C (G) is an involutory Hopf algebra. When G — Z 2 we obtain the toric code model. 
Thus the quantum double Hamiltonian realizes a representation of the quantum double of these 
involutory Hopf algebras. This is a general recipe to construct quasi triangular Hopf algebras 
from a given Hopf algebra. These quasi triangular structures are governed by a R-matrix which 
form a representation of the braid group in two dimensions and hence help us find anyonic 
solutions [21]. In [T7] we embed such quantum double models in a bigger parameter space 
namely the full parameter space of lattice gauge theories in three dimensions. For particular 
choices of the parameters we obtain the solvable models which describe the same phase as the 
QDM phase of Kitaev. For more general parameters we go away from this phase in a manner 
similar to the effect of perturbations which can be thought of as magnetic fields in the simple 
cases of the group algebra. By studying it in a bigger parameter space new topological and 
quasi-topological phases were also shown to exist H3 for these lattice gauge theories. These 
phases were understood by analyzing the Hamiltonians derived from the transfer matrices of 
these lattice gauge theories. These were constructed by carrying out Kuperberg’s prescription 
on a triangulated manifold with boundary. These transfer matrices are a product of local 
operators acting on the vertices, plaquettes and links. The vertex and plaquette operators 
are precisely the ones appearing in the QDM Hamiltonian |T8j [T9] . The ones acting on the 
links are like the perturbations to the QDM Hamiltonian. The trace of these transfer matrices 
coincide with the Kuperberg 3-manifold invariant. In this regard mi is a significant extension 
of Kuperberg’s constructions as it is not limited to just obtaining 3-manifold invariants. By 
also parametrizing the transfer matrix we go away from the topological limits and recover 
the topological invariants for special parameters which incidentally also give the topologically 
ordered models. 

The QDM and the string-net models are all examples of Hamiltonian realizations of lat¬ 
tice gauge theories based on involutory Hopf algebras and quantum groupoids or weak Hopf 
algebras respectively. The gauge fields in these examples live on the edges of the three 
dimensional lattice. These are examples of long-ranged entangled (LRE) phases which form a 
large chunk of known topologically ordered phases [23]. The “other” chunk of known topological 
phases are the short-ranged entangled ones. These include the Symmetry Protected Topolog¬ 
ical (SPT) phases which are interacting bosonic phases which have edge states protected by 
a global symmetry group m- Conventional topological insulators [25] are non-interacting 
fermionic versions of these phases. Such phases are known to exist in all dimensions, however 
the physically interesting ones are the ones in one, two and three dimensions. Exactly solv¬ 
able effective models describing such phases have been studied in [26, ,27]. Gauging the global 
symmetries of these systems lead to fractional excitations in the bulk. These help mimic the 
behavior of fractional quantum Hall (FQH) states [23] . 

Apart from these two broad classifications of topological phases, the LRE phases were further 
subdivided when they are decorated with an additional symmetry giving rise to symmetry 
enriched topological (SET) phases [2911301 131] . These phases were studied in the mathematical 
framework of G'-crossed braided fusion categories [32]. Exactly solvable lattice models for 
realizing SET phases with global symmetries were done in [33, 3T] [35]. Gauging these global 
symmetries leads to new topologically ordered phases starting from the parent topological 
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phase. Such phases were obtained locally in exactly solved Hamiltonians [?, ST* :3H] • They were 
realized in more realistic systems, such as bilayer FQH states, by introducing dislocations and 
the accompanying branch cuts [39]. 

A common feature of the above systems is the existence of a global symmetry. In the case of 
the SPT phases these global symmetries are realized by using matter fields living on the vertices 
with no link degrees of freedom [23] and for the SET phases the global symmetry is realized 
by placing the system on an enlarged lattice [33] (34- ;35] or by introducing dislocations in the 
lattice |39j. Thus it is desirable to go beyond the pure gauge systems in HD by implementing 
an additional global symmetry which may help one obtain the SET and SPT phases in a 
systematic manner. We do this by extending the construction in na by including matter fields 
on the vertices and construct the corresponding transfer matrices of systems with gauge and 
matter holds. By using this transfer matrix we obtain exactly solvable models of such systems 
which also includes the interaction between the gauge and matter fields. These Hamiltonians 
are frustration free HU and possess both a global symmetry, in the matter sector and a local 
symmetry in the gauge sector. They are frustration free in the sense that the Hamiltonian is a 
sum of terms such that the ground states of the full Hamiltonian are the lowest energy states 
of each individual term. The global ground states are also local ground states and hence there 
is no frustration or energy increase when all the terms in the Hamiltonian are included. A 
Hamiltonian which is a sum of commuting projectors satisfies these conditions and we shall see 
that the Hamiltonians constructed in this approach are of this form. 

Before we provide a description of the construction of these transfer matrices we briefly 
write down some properties of the systems obtained from these transfer matrices. As noted 
before the Hamiltonians are exactly solvable and the spectrum is gapped. The main input in 
building these systems is the gauge group and the matter module which is acted upon by the 
gauge group and the choice of the representation of the gauge group on the matter module. 
Depending on this choice we have can obtain a variety of systems. A common feature to all of 
them is the ground state degeneracy which is no longer a topological invariant like in the pure 
gauge case CD but rather develops another kind of degeneracy which depends on the number 
of equivalence classes of the matter module under the gauge action. Depending on the choice of 
the action of the gauge group the system can also develop a topological degeneracy in addition 
to the one coming from the different equivalence classes. These systems have excitations in both 
the gauge and the matter sectors. Unlike the deconhned excitations of the LRE phases, only 
the charge excitations are deconhned in this model. Depending on the choice of the action of 
the gauge fields on the matter degrees of freedom the system may or may not have deconhned 
huxes. When the huxes are confined, they are done so by string tension terms provided by 
operators acting on links as we shall see soon. However when these systems are placed on a 
manifold with a boundary the huxes get deconhned and become edge modes which are protected 
by an energy gap provided by the link operators. This feature will be explained with a specihc 
example in the text. This property of conhned bulk excitations and deconhned excitations on 
the edges has been observed in conhned Walker-Wang models mi elaborated in [421 SSI S3]. 
The conhned Walker-Wang models are systems in 3D, which when placed on manifolds with 
boundaries result in topologically ordered surface states. The examples in this paper achieve 
this in two dimensions through gapped edge modes. 

We briehy describe our method to construct the transfer matrix of a lattice system with 
gauge and matter helds in two and three dimensions thus corresponding to quantum lattice 
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models in one and two dimensions respectively. For convenience we only work with the case of 
finite groups though our methods can be extended to the case of continuous compact groups and 
in general to involutory Hopf algebras. Thus our focus is only on models which are pertinent 
for condensed matter systems. In particular we will be interested in finding exactly solvable 
Hamiltonians which act as effective descriptions of interacting systems of gauge and matter 
fields in one and two dimensions. We will explain the procedure for the three dimensional case. 
The procedure can be easily followed in the two dimensional case as well and is included as an 
appendix to this paper. 

The construction proceeds as follows. The data needed to define a partition function is a 
triangulated, closed 3-manifold which we take to be of the form E x S' 1 , where E is the two 
dimensional spatial slice, an involutory Hopf algebra A which make up the gauge degrees of 
freedom located on the edges of the triangulated 3-manifold and a n-dimensional vectorial space 
H n , which carries a representation of and has a co-structure and an inner product. The 
elements of H n make up the matter fields which are located on the vertices of the triangulated 
3-manifold. The invariant is built by contracting local weights associated to vertices, edges 
and faces. These local weights are built out of the structure constants of the involutory Hopf 
algebra A and the co-algebra H n . These structure constants from A are given by tensors m c ab , 
and S h a . They correspond respectively to the multiplication (product) map, coproduct map 
and the antipode. The structure constants from H n are given by tensors and /i“ ;3 where the 
former is similar to the coproduct for coalgebras and it arises from the costructure of H n and 
the latter is the map which says how the gauge fields act on the matter fields, that is it is the 
map which furnishes a representation of A on H n . In the structure constants listed the Latin 
alphabets index the gauge degrees of freedom and the Greek alphabets index the matter degrees 
of freedom. The scalar constructed out of these structure constants can be thought of as the 
result of contracting a three dimensional tensor network where we now consider the local weights 
as local tensors associated to the various parts of the three dimensional lattice. The partition 
function Z is parameterized by z and z*, elements of the centers of the algebra A and its dual, 
my an element of H n and G the inner product in H n . Thus we obtain Z ( A , H n , z, z*, my , G) 
which is in general a scalar and need not be topologically invariant. 

Obtaining an operator from this scalar is a natural step when we carry out the above 
procedure on a 3-manifold with boundary. This operator is precisely the transfer matrix as it’s 
trace is the partition function. The 3-manifold we consider is E x [0,1] where as before E is 
the two dimensional spatial slice and the unit interval is along the third direction which we can 
think of as the Euclidean time direction. Thinking of the partition function as a contracted 
tensor network we now have non-contracted indices on the spatial slices E x {0} and E x {1}. 
This results in the transfer matrix U for the lattice system. However since we distinguish spatial 
and “time” directions in this construction we have more parameters for U. This means that the 
weights associated to the two directions are not necessarily the same. Denoting the space and 
time directions by S and T respectively we now have the fully parametrized transfer matrix 
as U (A , H n , zs, zt, z$, my, Gs, Gt )• We show how this global operator can be written as a 

1 In other words, H n is a left H-module. 


5 



product of local operators acting on vertices, links and plaquettes. We obtain 

U ( zs , z T , z* s , Zy, m v , Gs, G T ) = n B r {z s ) n Q(G S ) n 

p i i 

I] Vv(G T )Q v (m v )l[A v (z* T ) (1) 

V V 

where v, l and p denote the vertices, links and plaquettes respectively. 

This operator is very general and encompasses a wide variety of interacting lattice models 
with gauge and matter fields. However for specific values of the parameters, namely 


zt = V , 



m v = Hn e , 


{G T )a/3 — , 

where p and A e are the unit and co-unit of the algebra A and Hn e is the counit of H n , we obtain 
exactly solvable models. The set of parameters we will work with in this paper are z$, z^ and 
Gs which give us 

u(z s ,4,Gs)= n^n^n^). ( 2 ) 

p l v 

These give us models which are similar to the toric code in the sense of a frustration-free 
Hamiltonian made up of commuting projectors and having long ranged entangled ground states, 
but including matter fields. This Hamiltonian is given by 

# = -5>”-E b »-E c ‘ ( 3 ) 

v p l 

where the vertex operators are gauge transformations just as in the QDM case, the plaquette 
operators measure gauge fluxes around a plaquette, and the link operator which is a gauged 
term with Potts-like nearest neighbor interactions. In other words, Ci couples the gauge con¬ 
figurations at the links with the Potts spins located at the vertices. 

The contents of this paper are organized as follows. In section 2 we describe the construction 
of the transfer matrix in 2 + 1 dimensions in a systematic manner starting from the method to 
construct the weights/local tensors using the structure constants of the input algebra for the 
gauge algebra and its corresponding vector space carrying it’s representation. The partition 
function resulting from this assignment is briefly sketched before the transfer matrix is obtained 
from this using the splitting procedure. The splitting procedure leads to the transfer matrix 
written as a product of local operators. All these local operators are obtained at the end of 
section 3 as tensor networks, written down in the Kuperberg notation ffl- The algebraic 
expressions for these operators are given for the gauge algebra being the group algebra and the 
vector spaces being the ones carrying their representations, in section 4. The input algebra can 
be more general than group algebras with group algebras being most relevant for condensed 
matter physics. Some examples of exactly solvable models obtained from this construction are 
presented in section 5. A brief outlook makes up section 6. There are two appendices intended 
as a supplement to the main text. Appendix A furnishes the details of the basic input for our 
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construction. In appendix B we show the same construction in 1 + 1 dimensions which produces 
quantum models in ID. These can be used to obtain spin chain models and hence completes 
our formalism for physically relevant instances. Finally in appendix C we look at two examples 
of exactly solvable models in one dimension and find them to have interesting properties which 
warrant further study. 


2 Constructing the Transfer Matrix 


In order to define our model we define a lattice £ which is a discretization of some (2 + 
1)—manifold £ x S 1 , where £ is some compact 2—manifold such that <9£ = 0. For the purpose 
of this work it is enough to consider £ a square lattice as the one shown in figure 1(a) This 


lattice is constructed by gluing vertices, links and faces as shown in figure 1(b). The model 


Z=7I 



(a) A piece of a (2 + 1)D square 
lattice. 



(b) Small pieces of which the 
(2 + 1 ) lattice is made of. 


Figure 1: Square lattice and its pieces. 

we build has degrees of freedom associated to gauge fields (living on the links) and degrees of 
freedom associated to matter (living on the vertices). They are quite general in the sense that it 
can support models such as the quantum double models which includes the toric code, besides 
models in other phases of matter. The Hamiltonian operator is the one which propagates the 
states along the time direction and it is made up of a set of projectors operators which act on 
specific parts of the lattice. In this paper we will show how to get the Hamiltonian operator 
in the language of tensor networks. The way we proceed is very similar to the way we have 
done so in dt We start with a partition function Z written in the formalism of the state sum 
model which leads to a one step evolution operator U such that H 

Z (^Sj z t, z£, z t> m v, Gs, G t ) = tr (U N ) , 

where Zg, Zt, Zg, are elements of the center of algebra and co-algebra which play role of 
parameters of the model, the subscript S and T means spacelike and timelike parameters while 
N is the number of steps in the time evolution. From this transfer matrix we shall be able to 
get a Hamiltonian H by taking its logarithm, i.e., 

U = U (zg, zt, Zg, z t, m v, Gs, Gt) — e At H . 

2 For brevity in the partition function expression we omit the dependence on the lattice (£), 
on the algebra (H) and on the H n space. So by Z (zg, Zt, +31 m v, Gs, Gt) we actually mean 
Z (£, A, H n , z s , z T , zs> z t, m v, G s , G T )- 
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In the next section we will define the function partition function Z (zg, z^, z^, z^,uiv,Gs,Gt) 
and relate it with a one step evolution operator. The partition function is made of a bunch of 
tensors associated to the vertex, links and faces of the lattice, all of them contracted with each 
other as we will see next. 


2.1 The Partition Function 


The partition function of the model is defined on an oriented square lattice £. We choose a 
square lattice for convenience. It can be easily defined on an arbitrarily triangulated lattice. 
The lattice is made up of vertices, links and faces glued together as shown in figure 1(b) A 


tensor is associated for each vertex, link and face of £ and contracted according to the gluing 
rules described below. In the next section we will show how the tensors we use to define the 
partition function are related with the structures constants of A and H n . 


2.1.1 Associating Tensors to the Lattice 

The procedure outlined helps define a function £ —>■ C. 


£ 


= 


— “1 -- “3 


(a) The tensor M aia2a3Q4 associated to a plaquette of 
the lattice. 




<*6 
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(b) The tensor T 011012013014015016 associated to a vertex of the 
lattice. 



(c) The tensor associated to a link of the lattice. 


Figure 2: The indices labelled by Latin letters (black arrows in the Kuperberg diagram) stand 
for the gauge fields while the indices labelled by Greek letters (green arrows in the Kuperberg 
diagram) stand for matter fields. 


A tensor M, 


a l a 2 a 3 a 4 


is associated to each face. The four indices label the four links glued to 


this face. This is shown in figure 2(a) The indices are ordered counter clockwise according to 


orientation of the plaquette. Each vertex is glued into six different links and hence the tensor 
associated to it is denoted T 01 " 203 " 4 " 5 " 6 as shown in figure 2(b). Finally each link is glued 











into four faces and two vertices, and so the tensor associated to a link is L abcda p as shown 
in figure 2(c) Note that the indices labelled by Latin letters (black arrows in the Kuperberg 


diagram) stand for the gauge fields while the indices labelled by Greek letters (green arrows 
in the Kuperberg diagram) stand for matter fields. In figure [2] the dotted lines mean that 
something is going to be glued to it. 


2.1.2 Contraction Rules 

To take care of the orientation of the lattice we still need two more tensors which will play 
the role of fixing orientation of the gauge part and the matter part. The one associated to the 
gauge part is the antipode map S of the Hopf algebra A while the one associated to the matter 
part is some bilinear map G : H n <g) H n —y C which we will describe later. The Kuperberg 
diagram for them are the ones shown in figure [31 

G a/3 = a 0 

(b) The orientation tensor 
associated to the matter 
part. 

Figure 3: The tensors which take into account the orientation of the lattice. 


S b a = 

(a) The orientation ten¬ 
sor associated to the gauge 
part. 


The contraction tensors with the contraction rules are shown in figure [H For the gauge 
part (contraction of a plaquette with a link) we use the antipode tensor when the orientation 
does not match or we contract the tensors M aia2a3a4 with L abcda p directly as shown in figure 


' abcdoL 


4(a) For the matter part we contract the green arrow which is going out of the tensor L a -^ 

using the tensor G a p while the green arrow coming in the tensor L abcda p we contract directly, 


as shown in figure 4(b) 





(a) The contraction rule for the tensors in the gauge sector. 


(b) The contraction rule 
for the tensors in the mat¬ 
ter sector. 


Figure 4: Contraction Rules. 

The partition function is defined as the contraction of all these tensors associated to the 
vertices, links and faces. Since the lattice we are considering has no boundary, dC = 0, there 
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will be no free indices remaining after the gluing resulting in a scalar. As we will see in the 
next section the tensors M aia2a3a4 and L abcda p depend on the structure constants of A, H n , 
elements of the center of the algebra A and its dual, elements of H n . With these parameters 
the partition function can be written as 

Z (£, z s , z T , zl 4, m v , Gs, G T ) = ]} J] Sj f] G af , 

p v loo 

where the products runs over the orientations (o), the plaquettes (p), links (/) and vertices (v) 

of C. 


2.2 The Transfer Matrix U from the Partition Function Z 

Having defined the partition function we can now define a one step evolution operator U such 
that 

Z = tr (U N ) . 

The way to do this is by looking at the lattice in which the partition function is defined on 
and to take slices of it in the time direction. We now introduce a new point of view where 
we look at the lattice as a contraction of a bunch of tensors as defined above. As we will see 
sometimes it will be convenient to think of the lattice as a tensor network and sometimes more 
convenient to think of it as a gluing of vertices, links and faces. In figure 0 the black and 


Z {£,■•■) 





Figure 5: The partition function as the trace of a product of a one step evolution operator. 


green arrows represent the free indices (gauge and matter degrees of freedom) of the one step 
evolution operator. The one step evolution operator is the one shown in figure [6J, it is made of 
a bunch of boxes without their caps glued one besides the other. 
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Figure 6: The one step evolution operator U. 


3 The Transfer Matrix as a Product of Local Operators 

We will see that the one step evolution defined in figure |H] is a product of local operators. 
These local operators act on vertices, links and plaquettes of the lattice, and hence are called 
vertex, link and plaquette operators, respectively. But before we split it we must have a better 
definition of the tensors used to build it. In other words we need to know how they are related 
with the structure constants of A and H n . In the next section we will define their algebraic 
structure which is a very important part of the model. 


3.1 The Structure Constants 

Let A be an involutory Hopf algebra with a product m : A® A —> A, a co-product 5 : A —* A® A 
and a involutory map S : A —> A. Its structure constants are the tensors shown in figure [7] 


—*- 





(a) The 

(b) The 

(c) The an¬ 

multipli¬ 

co-product 

tipode map. 

cation 

map. 


map. 




Figure 7: Structure constants of the involutory Hopf algebra A. 


Let H n be a module of A which has a semi-simple co-algebra structure t : H n —> H n <g) H n 


We define the action of A on H n by the map /i : A ■ 


H n —> H n in the following way 


/i : a <g) v •—> /i(a) > v . 


The structure constants involving H n are the ones shown in figure 8(a) and 8(b). We also show 
the action of the algebra A on H n with respect to the product of the algebra A, in other words 
we want this action to be a homomorphism, i. e., fi(ab) > v = fi(a) > (/i(6) > v) Wv as shown 


in figure 8(c 
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(a) The action map of A (b) The co- (c) The algebra action on the module is a 

on H n . structure homomorphism. 

map in H n . 


Figure 8: Structure constants involving the H n space. 

3.2 Building the Tensors 

The first tensor we will define is the tensor M aia2a3a4 . Its definition is shown in figure El Due to 
the associativity of the algebra the tensor M aia2a3a4 is invariant under cyclic permutation of its 
indices and can be written in different ways as shown in figure EJ The element z is an element 
of the center of the algebra. 



«2 


Figure 9: The definition of the tensor M aia2a3a4 . 


The tensor T aia2a3a4asa6 is made of the co-structure map in H n as shown in figure [lOj As 
the co-product is co-commutative this tensor is invariant under interchange of a pair of its 
indices. Like the tensor M aia2a3a4 this tensor can also be written in different equivalent forms 
as shown in figure [TUI 



Figure 10: The definition of the tensor j 101101201301 * 015016 . 


Finally the tensor associated to the links given by L abcda p is defined in figure fill As we can 
see it is made up of the algebra action on H n given by the /i map and the tensor A aia2a3 “ 4a5 
which is also defined in figure [IT] The element is an element of the co-center of the algebra 

A*. 
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0,4 


Figure 11: The definition of the tensor L abcda p. 


3.3 Splitting the One Step Evolution Operator 

In order to make the procedure clear we will represent the spacelike part as being made of 
spacelike plaquettes and spacelike links and the timelike part as being made of timelike plaque- 
ttes and timelike links. The splitting will be done in two steps: first we will split the tensors 
associated to spacelike links and vertices and then we split the tensor associated to the timelike 
plaquettes. These two steps show how the one step evolution operator can be written as a 
product of local operators. Although the model does not depend on the orientation of the 
lattice we have to set up some orientation in order to use the contraction rules described above. 
So we consider the orientation shown in figure fl2l 



Figure 12: The orientation of the lattice considered. 


3.3.1 Splitting the Spacelike Part of the Tensor Network 

Consider a spacelike link which is shared by two adjacent plaquettes and two vertices as shown 
in figure [[3j 

The tensor network associated to the figure [13] is shown in figure [14l In the picture tip and 
til stand for timelike plaquette and timelike link respectively. 

As seen in figure [TT| the tensor L abcda p can be splitted in terms of the structure constants 
representing the co-product of the algebra. Thus the diagram in figure dH reduces to the one 
shown in figure [15] 

To make things clear we write down the same tensor network of figure [15] in figure [I6l where 
we have separated the tensors associated to each one of the spacelike plaquettes. Each leg of 
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Figure 13: A spacelike link shared by two adjacent plaquettes and two vertices. 


out 



Figure 14: Tensor network associated to the picture in figure fl3l 


out 



Figure 15: Tensor network associated to the picture in figure [13] with the tensor L abcda p splitted 


the tensor M aia2a3a4 will be contracted to a tensor 8^ d (directly or indirectly by the antipode 
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map) and that will lead to the plaquette operator as we shall soon see. 


out 



(til) 


Figure 16: The tensors network in figure fl5l rewritten for clarity 


We now split the tensor associated to each vertex in figure [[6j For that we just use the 
definition of the tensor T 0101201301 * 015016 in figure [TOl The splitting shown in figure [T7] has to be 
done for the all the vertices of the lattice, but here for simplicity we are illustrating just for one 
single vertex. 



Figure 17: Splitting the tensors associated to a vertex. 


The tensor network in figure [16] is rewritten in a separated way as shown in figure [TT1 After 
all this splitting we see that the spacelike tensor network part can be written as a product of 


two blocks of tensor networks which are the ones shown in figure [TH] The first one (figure 18(a) 


is the operator called the plaquette operator. The second one in figure 18(c) is called Ki and 
it will be contracted with the timelike tensor network part to build the others operators which 
make up the transfer matrix. 
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©—*■ <5 


B p 




~L —- 


(tip) 


T 


(a) The plaquette operator as a ten¬ 
sor network. 


(b) The 
plaquette p 
where the 
operator B p 
acts on. 


(c) The tensor network 
which will be contract with 
the timelike part of the 
remaining tensor network. 


Figure 18: The spacelike part of the tensor network can be written as a product of these two 
blocks of tensor networks. 


3.3.2 Splitting the Timelike Part of the Tensor Network 


We now split the timclikc part of the tensor network. For that consider a timelike plaquette 
which is shared by two timelike links and one spacclike link, as shown in the figure 19(a) The 


spacelike link on the top is the one which has already been splitted and it led to the tensor 


network attached to the timelike plaquette as shown in figure 19(b) 



in 


(a) A timelike plaquette of the transfer ma¬ 
trix. 



(b) A timelike plaquette with the tensor 
network associated to the spacelike link at¬ 
tached above. 


Figure 19: Timelike part of the transfer matrix. 
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The tensor network associated to the timelike plaquette contracted with the two timelike 


links is the one one shown in figure 20(a) and it can be easily changed to the one in figure 20(b) 
by using the definition of the tensor M aia2a3aA in figure [9j 


I 

S 

\ 

-S — 


I 


(a) The tensor network associ¬ 
ated to the timclikc part of the 
diagram. 




-m- 


m 


S 

\ 

M 

t 

► m - 


t / 


i 


m m 

(b) The tensor M aia2a3ai splitted in 
terms of the structure constants of the 
algebra. 


Figure 20: Splitting the tensor network of the timclikc part of the transfer matrix. 


As before we rewrite the tensor network in figure 20(b) as the one in figure I2T1 



Figure 21: Rewriting the tensor network in figure 20(b) with the tensors separated. 


In order to keep it clear in our mind where are these tensors acting on one should take a 
look at the figure 19(b). Note that each and every timclikc tensor L abcda p will have a tensor 


m c ab contract with each of its gauge legs (directly or indirectly with the antipode). It lead to 
an operator A v for each timelike link of the lattice, or in other words, for each vertex spacelike 
vertex v. This operator is the one drawn in figure 22(a). The tensor in highlight in figure 


is now attached to the tensor network on the top in figure 19(b) and it give us a new operator 


Ki which acts on the links, see figure 22(c) 


The operator Ki is the same on which appears in the (1+1) dimensional models. In the 
appendix we show how the transfer matrices of such a models can be obtained e also that the 
Ki operator can be written in terms of simplers operators. The transfer matrix we started with 
can now be written as 


U (zs> z t, 4 z t, m v , Gs, Gt) — B p (zs) A) A v (z^) . 

p l V 

But the Ki operator, as can be seen in the appendix, can also be written as 

=n^(Gs)n( T '(^)AA(^))n^( G r)<9«K), 

III V 
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(a) The vertex operator. (b) The ver- (c) The tensor network of the link oper- 

tex v where ator. 

the vertex op¬ 
erator A v acts. 


Figure 22: Vertex operator and the operator Kj which acts on a link. 


where the operators Q v (m v ), Ti(zg) and Li(zt) will be explained in the next section. The 
operators Di and V v (Gt ) are explained in the appendix. The operator Ci(Gs ) is called the 
link operator and acts on a link and on the its vertex. This operator is defined below in figure 
23(a)] 


(— 



(a) The link operator. (b) The link l where the link 

operator Ci acts. 


Figure 23: The link operator. 


4 The Transfer Matrix U in its Final Form 

So far we have used a diagrammatic language where we associated tensors to the vertices, links 
and plaquettes of a triangulated three manifold to build a tensor network and showed that this 
results in a partition function for a closed manifold and a transfer matrix for a manifold with 
boundaries. We then used this transfer matrix represented by a tensor network to find out the 
local operators which can be used to piece up such a network. In this section we will write 
down the algebraic expressions for these operators and study their properties. For simplicity 
we choose the gauge algebra to belong to groups and so the states on the links are elements of 
the group algebra of a gauge group denoted by C(G). The matter degrees of freedom belong 
to a module of C(G) which we denote by H n , in other words H n is a vector space carrying 
an action of C(G). The formalism developed so far can be applied to any involutory G*-Hopf 
algebra and its corresponding module. 

The operators that make up the transfer matrix can be divided into the ones which act only 
on the gauge fields, the ones that only act on the matter fields and the ones that involve the 
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coupling between the gauge and the matter helds through the /i map defined in figure (TTTi) . We 
describe each set separately before we write down the full transfer matrix. 

In order to write an algebraic expression for the operators we derived in previous section we 
need to define three operators which are made of the structure constants of the algebra and the 
module, as shown in figure [2U These operators are linear on its parameters, in other words, 


L(z) = 

z 

1 

— *m —*- 

— ► 

R(Z) EE | 

z 

2* 

T(z*) = \ 

Q(V) EE J 

(a) left 
tion. 

multiplica- 

(b) right multiplica¬ 
tion. 

(c) gauge projector. 

(d) matter projec¬ 
tor. 


Figure 24: These operators are the ones which we use to build the vertex, plaquette and link 
operators. 


L(z) = Yl g z 9 L(4> g ) (the same for R(z), T(z*) and Q(v)) and they act on the vector basis as 

l h ) = 1 9h) , 

R{4>g)\h) = | hg) , 

T(* g )\h) = S(g, h)\h) , 

Q(Xi)\j) = s{i,j)\j) . (4) 

where we have used the bra-ket notation for elements of the basis (|g) := (j) g and |i) = %*). 
Sometimes we use the short notation L 9 = L((f> g ), R 9 = R(4> g ), T 9 = T('F £/ ) and Q l = Q{xi)- 
Now using this operators in the next section we write the plaquette, vertex and link operators 
derived before in terms of them. There are three operators acting on the gauge sector alone. 
These are given by the plaquette operator and two operators acting on the qudits on the 
links. They can be deduced from their respective tensor network representation from the 
previous sections. We write down each of these in what follows. We will also see how the 
different parameters parametrizing the transfer matrix U (zs, Zt, Zg, z^,my,Gs, Gt) arise in 
the definition of these operators. This will also show how the transfer matrix depends on these 
parameters. 


4.1 The Gauge Sector 


The plaquette operator can be written down by looking at its tensor network representation 
given in figure (18(a)). We find 


b„ = Y. B l r (s) 

Ce[G] 

where C labels a particular conjugacy class from the set of conjugacy classes in G denoted by 
[G]. Each of is given by 

Bp = M E E *(01320304 9 ,1 oK" 0 Tff <8 T» 8 Tf;' (6) 

9^C {gi} 

where 1q is the identity of the group G. 
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While writing down the plaquette operator in Eq.(15j) we have made use of the fact that the 
tensor associated to the spacelike plaquette M aiQ2a3a4 contains a central element of C(G) as can 
also be seen in figure ([9]). For the plaquette operator in Eq.Q this central element zs, where 
S denotes spacelike, is given by 

z s = Y P c Y 4 ( ? ) 

ce[G] <?£C 

with the set {4> g } forming a basis of C(G).The action of the plaquette operator on the oriented 


square lattice is shown in figure 18(b) 

The operators form a complete basis of orthogonal projectors since 

BSB? = S(C, C')Bp , 


and also 




c 


One of the link operators acting on the qudit on the links is given by 

Ti(z's) = EW 

g&G 


( 8 ) 


where b g are constants and Tf is gauge projector defined in figure 24(c). These operators arise 
by the parameter on spacelike links given by z$ which is a central element of the dual of C(G). 
They can also be seen in figure (Hill . This central element is given by 

4 = Y b 9^9 (9) 

g&G 

with the set {r/^} forming the basis of the dual of C(G). 

The second link operator acting in the gauge sector is given by 

L j( z r) = Y a R Lj(z R ) (10) 

.ReIRR’s of G 

where a R s are constants and z R is given by 

1 


zr = 


|G*| 


Y Xr 4)4- 


(ii) 


g&G 


Xr4) is the character of the element g in the IRR R and |G#| is the number of elements with 
non-zero trace in the IRR R. The parameters z R is given by 

zt = Y^ arZr (12) 

R 

We could as well have defined the operator Rj(zt), but since z R is an element of the center of 
the algebra these two operators are the same, namely Lj(zt) = Rj(zt)- 

The operator X 3 in Eq. flTOl) is obtained through the parameter z R in the timelike plaquette. 
This is an element of the center of C (G) given by 

Z T = Y a R Y x 49)(t>g- (13) 

R g€G 

These operators exhaust the operators acting only in the gauge sector. We now turn to 
those acting only in the matter sector. 
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4.2 The Matter Sector 


There are two operators which act only in the matter sector. One of them is obtained by 
parametrizing the vertices by my, an element of the module H n and the other by parametrizing 
the timelike links with the inner product Gt- The operator parametrized by my is given by 

n 

Q v (m v ) = ^TciQ l v (14) 

2 = 1 

where c t are constants and Q l v is given by 

Qv\j) = koti) (15) 

where Xi,j are elements of the module H n . For the definition of Q v in Eq. (fT4j) the element my 
is given by 

n 

my = ^2,dXi. (16) 

2=1 

As the module only has a co-structure this operator is very similar to the coproduct map 
in the gauge sector. In fact for the module H n with symmetry group Z n we can identify Q J V 
with the T uJ operator of the gauge sector for C(Z n ) with u = e4r and j G (0, • • • ,n — 1). If 
we label the elements of the group Z n as {c o k ] k G (0, • • • , n — 1)} we can then write Q l v as 

1 n— 1 

= (17) 

1 k =0 

where is the character of the element uj k in the IRR of Z n labelled by c o l . There are n 

such expressions corresponding to the n IRR’s of Z n . 

The operator Z k in Eq. flTTll is a generator of Z n and is defined by 

Z v \u k ) =u k \oj k ). (18) 

We can also consider modules H n with other symmetry groups, especially non-Abelian 
groups, and in this case the expression for Q l v is different from the one given by Eq. flTTlh The 
symmetry groups of the module H n can also be thought of as global symmetry groups of the 
system. 

We now consider the operators acting on both the gauge and matter sectors. 


4.3 The Gauge + Matter Sector 

There are two operators acting on both the gauge and matter sector. They are the operators 
coupling the two sectors. These operators are the vertex and link operators given in figures 
(22(a)) and (23(a)) respectively. 

The vertex operator A v is given by 


Ay = a gAy 
g&G 


( 19 ) 
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where 


Al = /i„(0 s ) ® L‘ h ® if, ® R %'' ® R%' (20) 

with n v (4>g) being the representation of the gauge group on the matter field located on the 
vertex v. The single qudit operators Lf and Rf, acting on the gauge fields located on the links, 
are given by Eq. (1Tj) . 

The vertex operator A v is obtained by using the parameter z £ in the transfer matrix given 
by 

4 = (2i) 

g&G 

This is a parameter living on the timelike links of the 2 + 1 dimensional manifold. 

For particular choices of the parameter z £ we obtain the set of orthogonal vertex operators 
which are projectors. These parameters are given by 


z 


* _ 

T — 


l 


J2xr(9)^ 9 

g&G 


( 22 ) 


for the different R 1 s in the set of IRR’s of G. They result in the following set of vertex operators 


A 


R 

V 


1 

i^f 


^2xr( 9)A 9 v . 

g£G 


(23) 


This vertex operator is very similar to the one defined for the quantum double models m 
HEIDI with the addition of the representation of the gauge group acting on the vertex part. 
Thus it can still be thought of as a gauge transformation just as in the case of the quantum 
double models. 

The other operator acting on both the gauge and the matter sector is the link operator 


represented as a tensor network in figure (23(a)). In terms of operators this operator is given 
by 

c,= V G s {p(4>,)Xv„X,Ml: l Tt‘Qll' (24) 

Xvi > Xv2 1 + 


where the matrix Gs implements the inner product between the vectors in the module H n . 

This inner product G can also be thought of to be represented by the following tensor 
network shown in figure (1251) . where T~ l is a tensor such that ( T~ l ) a pT l3l = < 5 (ck, 7 ), in other 


G- - = - -R - -T- 1 - - R 


Figure 25: The tensor network representation for the inner product G in terms of the matrix 
R. 


words, (T “ 1 ) a /3 = S(a,/3). 

The matrices R in this definition can be thought of as those performing basis transformations 
on the elements of the module H n . If the matrix G is chosen to be identity, it is the same as 
choosing R to be identity. For other choices of R we can find the operators orthogonal to C). 
This will soon become clearer when we illustrate using examples. 
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We thus have the operators forming the transfer matrix of a lattice theory with gauge and 
matter fields. To sum up these operators include the vertex, plaquette and link operators along 
with the single qudit operators acting in the gauge sector given by Eq.flH]) and Eq. flTO and 
in the matter sector given by Eq. 01 4 1) . We now proceed to study the algebra between these 
operators. 

4.4 Algebra of the Operators 

The remaining basic operators in the theory are given by Lf, Rf, Tf and Q l v . Clearly Tf and Q l v 
are projectors and Q l v commutes with the remaining operators as they act on different sectors. 
We have the following relations 


LfLf 

_ L?(92 

(25) 

D9l d92 
K i K l 

_ ft92gi 

(26) 

T ,, L ,, 

_ Jj}lrjl9\ 92 

(27) 

rjl92 J^gi 

_ j^9i r j'929i 

(28) 


Using these relations it is easy to see that the vertex and plaquette operators given by Eq. (TT9j) 
and Eq. (EJ) commute. This computation is exactly similar to the computation used to prove 
the commutativity of these two operators in the case of the quantum double models. This is 
true because the plaquette operators in our case of theories with gauge and matter fields is left 
unchanged with respect to the case with pure gauge fields. 

The link operators C) trivially commute with the plaquette operators as both are diagonal. 
This can be seen from their respective definitions given by Eq.(J2U) and Eq.(f>]). The only thing 
that we have to prove is the commutation between the vertex operator and the link operator 
in their common region of support. 

Consider the action of the operators A v and Q on the gauge and matter degrees of freedom. 
The proof then goes as follows for their action on the common support, we distinguish two 
cases as the vertex operator action depends on the orientation of the lattice, the first case being 
the one when the vertex operator acts on the vertex at the left of the edge l, this is 

-'A 1 Cj |. • • , Xvi ) (fill Xv21 • • • } {l^(,^fil)Xv 17 Xv2^) G- 4^1 I" ■ ■ ) Xvil ( fill Xv21 • • • ) 

{/^{(fifijXv 1 1 Xv2 )g I ■ • • j l^{(fig)Xv 1 1 (filtfig- 1 1 Xv21 • • • ) i ( 29 ) 

whereas 

C'l-Avi I • ■ ■ i Xvi i (fill Xv2 !•••'} f-U | ■ • ■ i /-((.(fig^Xv 1 1 (fil&g— 1 1 Xv2 i • • • ) 

= (/^{(filtfig- 1 ) f Jj ( < fig)Xv i) Xv 2 )g |- • • j l^{(fig)Xvii (fil^fig- 1 ! Xv 2 i • ■ ■) i (30) 

since the module map is a group homomorphism we note that: 

tl{(fil(fig-l)v((fig) = fl((fil)fl(^>g-l)fl((fig) 

= M)K(fig-i(fig) 

= M0A 
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thus, the coefficient on the right hand side of eq. (!291) now is the same as the one in eq. (j3U|) . 
hence [A^, C)] = 0. Let us now consider the case when the vertex operator is acting at the 
right of the edge Z, i.e. 

^V2 I ■ ■ ■ 7 Xv\ 7 0/ 7 Xv2 1 ■ • ■ ) (M (0l)Xvi 7 Xv2 ) gA V2 | • • • , Xvi 7 4*17 Xv2 7 • • • ) 

(a(00Xi; 15 Xv 2 )g | ■ • • 7 Xvi 7 4 ) g4 > h M(4 > g)Xv2 7 ■■■) 7 (31) 

while 

2 I ■ ■ ■ 7 Xv\ 7 0/ 7 Xv2 7 ••■ ) 10 | ' ' • 7 Xvi 7 4 > 9 4^l 7 1^(4*<j)Xv2 7 • • ■ ) 

{f^{4 ) g4 > l)Xv i ) f Ji ( 4 > g')Xv 2 ) G \ ■ ■ ■ 7 Xvi 7 4 ) g4 ) h A(0 9 )W 2 7 ■■•) 7 (32) 

for the commutation to hold we require: 

(M0 9 )Xa, Xp)G = (Xa, ^{4>g) ] Xp)G7 (33) 

this means the representation map /i is unitary, equivalently 

M0 9 ) f = ^(05-0- (34) 

therefore the right hand side of eq. fEEZ]) reduces to that of eq. (TBTjl . Thus, [A -\ yi , C)] = 0. So, we 
just showed that [A 9 vl Ci] = 0 for any 0 g e C(G) and for any vertex v of the lattice, thus it 
follows that [A v ,Ci\ = 0. 

Moreover we only consider /As which permute the basis elements of the matter module. 
Under this condition (fi(4>i)-Xvi \Xv 2 ) is real which makes the link operator Ci hermitian. This 
also ensures that the vertex operator is hermitian as the gauge part of the vertex operator is the 
same as the one in the QDM of Kitaev. The plaquette operator is the same as those appearing 
in the QDM of Kitaev. Thus all the operators in the Hamiltonian are hermitian making the 
evolution unitary. 

4.5 The Transfer Matrix 

We are now in a position to write down the full transfer matrix of a lattice theory with gauge 
and matter fields as a product of these local operators. Thus we have 

U{z s ,z T ,z* s ,z* Tl m v ,G S 7 G T ) = B p (z s ) J] G(G 5 ) ]J (104) D^zt)) 

p i i 

JJ V v (G T )Q v (m v )A v (z^) 

V 

The other operators namely B p (zs), Ti(zg), Q v (m v ), Li(zt), A v (z £) are given by Eq. (J5]) , Eq.flH]), 
Eq. Oldp . Eq. fllOp . Eq. 01 9 p respectively. Through these operators the transfer matrix U obtains 
its dependence on zs, z* s , m v , zt and z respectively. 

The operator Ci is defined as a linear combination of orthogonal projectors. We write it as 

Q = 7Q + 7±C0 (35) 

where Ci is given by Eq. 11241) and C l 0 is the operator orthogonal to C0 The form of the 
orthogonal operator C l 0 depends on the module H n considered. In general there are several 
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operators orthogonal to the operator C). For the modules with symmetry group as Z n we can 
write down a general form for these orthogonal operators. If we denote the inner product in 
the module H n as (Xi\Xj) then we can write down the operator Ci as 


C'l I Xvi ) 0/ j XV 2 ) (/^(0z) -Xvi Ix^) \Xvi ) XV 2 ) ■ (36) 

We can now write n — 1 other orthogonal operators as 

C^XviAhXv*) = (M)-Xv 1 \X l \xv 2 )\Xv 1 ,(l)i,Xv 2 )-, *e{V--,n-l} (37) 

where X is the shift operator generating Z n given by 

X k \ Xi ) = | Xi+k). (38) 

The transfer matrix gets its dependence on G through the link operator Q. The inner product 
given by G in turn depends on the matrix R as shown by it’s tensor network representation in 
figure (1251) . 

Thus for each i in Eq. (1H71) we choose a different R matrix in the inner product G. For a given 
i in Eq. (l37[) R is given by (X*) 2 . This matrix R can be thought of as a basis transformation 
in the module H n . 

Having written down the general expression for the transfer matrix U of a lattice theory with 
gauge and matter fields we will now consider specific examples which illustrate the formalism 
developed so far. 

5 Examples 

We will consider two examples of Hamiltonians having long-ranged entangled ground states 
derived from the construction we have illustrated. The first example is the simplest Hamiltonian 
with the gauge group being Z 2 and the vector space carrying its representation being the two 
dimensional vector space H 2 . We will denote this model as H 2 /’L 2 . The second example is 
a slight modification with the vector space H 2 replaced by H 3 , a three dimensional vector 
space carrying the representation of the gauge group Z 2 . We denote this model as H 3 / Z 2 . We 
consider each of them separately in what follows. 

The ground state degeneracies in each of these cases is a topological invariant and can be 
computed numerically. There are no winding operators as in the toric code for the 77 2 /Z 2 case 
but we do have a winding operator for the H 3 / Z 2 case. We will write down the ground states 
in both the cases and give their tensor network representations. 

Both these examples have confined excitations apart from one deconhned excitation. We 
will briefly look at the vertex and gauge excitations in the H 2 / Z 2 case. The vertex excitations 
of H 3 / Z 2 were studied in detail in [35|. Here we will just look at it’s deconhned excitation. We 
then present a third example, the H 2 / Z 4 model where we use an action of Z 4 on H 2 such that 
one of the fluxes become deconhned even in the presence of the link operator Ci. 

The effect of confinement is due to the link operator Ci. A systematic way to deconhne all 
the huxes is by working with Hamiltonians which do not have the Ci operators but just the 
vertex operator A v and the plaquette operator B p . We then see that we have all the deconhned 
excitations of the quantum double models along with the vertex excitations due to the presence 
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of the matter fields. These models will comprise the fourth set of examples. This is the simplest 
and most straightforward way to deconfine all the fluxes in contrast with the H2/Z4 case where 
a clever choice of input data and structure constant helped us deconfine the fluxes. 

Finally we will show how to recover the quantum double Hamiltonians with only gauge 
fields by “switching off” the matter fields. These comprise the fifth set of examples. All the 
examples are constructed on the torus. 

Before we write down the models we make a few comments about the ground state degener¬ 
acy of these models. The models we will construct are obtained from the transfer matrix which 
has the following form 

^n^n^iiA (39) 

p l V 

which makes this operator a projector. Using the arguments of El we can show that the trace 
of this operator gives the ground state degeneracy. With a little more work it is possible to 
show that this is also a topological invariant. This requires the introduction of Pachner moves 
in the presence of matter fields. We reserve this discussion for another paper. 

The trace of the transfer matrix can be computed numericalljH and these results are used 
for the two examples to be discussed below. 


5.1 H 2 / Z 2 : 

The basis elements of the group algebra of Z 2 are denoted by the elements of the set 
and the basis elements of the module H 2 are denoted by {xi, X-i}- The element </>i acts as the 
identity element on the module H 2 . The element </>_ 1 on the other hand flips between the two 
elements of H 2 . That is the representation map /i(0-i) is given by 


M0-i)Ai = X-i (40) 

= XT- ( 41 ) 


As a matrix is the cr x Pauli matrix and /i( 0 i) is the two by two identity matrix. Hence¬ 

forth we will use the Pauli matrices in place of the representation matrices. 

The two orthogonal vertex operators are given by 


I 0 I 0 E 0 I 0 I + a x 0 af 1 0 af 2 0 erf 3 <8 erf 4 

2 

I 0 I 0 I 0 I 0 I — erf 0 erf 0 erf < 8 > af 3 0 of 4 


The plaquette operators are the same as in the toric code case and are given by 

I 0 I (g) I + erf ® erf 0 erf 0 erf 

- = -2- 

10 10 10 I - erf 0 crf 2 0 erf 3 0 crf 4 


(42) 

(43) 


(44) 

(45) 


As can be seen from these expressions the plaquette operator acts trivially on the matter sector. 
3 We thank Kazuo Teramoto for these computations. 
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V\ 

/ 

V2 

1 

1 

1 

-1 

-1 

1 

-1 

1 

-1 

1 

-1 

-1 


Table 1: Configurations of Xvn 4>i and Xv 2 which make the inner product in Eq. ffTO]) non-zero. 


Vi 

l 

V 2 

-1 

1 

1 

1 

-1 

1 

1 

1 

-1 

-1 

-1 

-1 


Table 2: Configurations of Xv 1 , and Xv 2 which make the inner product in Eq. (l47|) non-zero. 


There are two orthogonal link operators in this case. Their action on the basis states of the 
system are given by 


f-'Z \Xv\ ) 0Z) Xv 2 ) 
C'l \Xv\i4 > hXv 2 ) 


(h( 0 i) -Xv\ | Xv 2 ) Ikhi ■> fill Xv 2 ) 

-Xv 1 \& IXU2) |Xi>i ■> (frh Xv 2 ) 


(46) 

(47) 


We now obtain the matrix representations of these two orthogonal operators. This is done 
by summing the projectors onto all the configurations that satisfy the constraint given by the 
inner product. The set of configurations Xvn <fti an d Xv 2 which make the inner product for Cj 
non-zero is shown in table (P). The projector to these configurations is written as the following 
sum 


This reduces to 



(48) 

(49) 


In a similar way we can compute C t 1 which is the orthogonal projector to (7/. The set of 
configurations an d Xv 2 which make the inner product for Cf 1 non-zero is shown in table 

©. 
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The projector to these configurations is written as the following sum 


CT 1 = 


1 + a. 


VI 


1 - < 


l + af 


1 — a. 


V2 


a: 


VI 


1 + a. 


VI 


err 


l + af 


1 ~° Z i 


2 

1 ~ <2 

2 

1 + a V2 


1 + a. 


V2 


This reduces to 


c -i = 1 < 8 > 1 <8 1 - < 


err Q9 a. 


V 2 


(50) 


(51) 


Using the matrix representations for all the operators it is easy to see that they indeed 
commute with each other and they are all projectors. The model posses a global symmetry 
given by the operator ]f[ j; af where the product runs over all the vertices on the lattice. This 
can also be thought of as the symmetry group of the module H 2 . This comparison will be 
crucial later on. 


l + af 


l + af 


The operator Xi is given by 

Xi = oi 

The operator Z\ is given by 

Zi = b i 

The operator Q v is given by 

Qv Ci 

The connector operator V v is given by 

T/ W < 1 + (JX 

V v = cti 


+ CL —1 


1 -af 


(52) 


+ 6 -i 


1 — af 


(53) 


1 + af 


+ c_i 


1 - af 


(54) 


+ d —i 


1 — a? 


(55) 


The other connector operator Li is proportional to identity nzt Thus we have the full 
transfer matrix for the case. We can now use the definition of the transfer matrix 


U (z s , zt, z* s , z ti m v, Gs, G t ) = e H 

to find several Hamiltonians having H 2 and 7L 2 degrees of freedom. One such example which is 
exactly solvable is given by the Hamiltonian 


(56) 























We now study the ground states of this Hamiltonian. The Hamiltonian in Eq. (16911 is exactly 
solvable as the operators making up the Hamiltonian are commuting projectors. Since the 
eigenvalues of these projectors are 0 or 1 the entire spectrum of the Hamiltonian is known. The 
condition for the ground states in particular is given by 

A v \gr) = B p \gr) = C t \gr) = \gr ); \/v,p,l. (57) 

This gives the lowest energy eigenvalue as — (N v + N p + iVj), where N v , N p and iVj are the 
number of vertices, plaquettes and links respectively. Numerical computations of the trace of 
the transfer matrix in this case shows that the ground state degeneracy is one in this case. 
One such state satisfying the conditions of Eq. (157|) is given by 

I i>pi) = U B p U C t ® 1 |A S ) <8)„ |A m ) (58) 

p i 

where \ g = + 4>-i and X m = xi + X-i- The representation of this state as a tensor network 

is shown in figure (|26|) . 


I^> 



Figure 26: The tensor network representation of the state | ip p i) given in Eq. (l58|) . 


We can write another ground state for this model as 

|^)= n ^®*l^)®-IXx>. (59) 

V 

The tensor network representation of this state is shown in figure (J 27]). 



Ixi) % lxi> 

Figure 27: The tensor network representation of | ip v ) given in Eq. (I59[l . 


The two ground states | ip p i) and \i]) v ) are similar to the two ground states in the case of the 
toric code as discussed in HU- They can be thought as being written in the basis of a x and a z 
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respectively. From our numerical considerations we presume that these two states are the same. 
Another way to see this is that there is only one equivalence class under the gauge action, that 
is the representation vector space has a single orbit under Z 2 action given by H 2 itself. Another 
reason why we expect these states to be the same is that there are no winding operators in this 
case and we know from our argument in HU that using the presence of winding operators we 
can relate the two basis sets of states. 

In this model we have a third possibility for a ground state, formed out of a mixture of 
states from the a z basis and the dual basis, given by 

\i/>vi) = n^Il^ l^ 1 ) 1^™)- (60) 

V l 

As the ground state degeneracy is one and there are no winding operators we again expect this 
state to be the same as \i]) v ) and \ip p i). 



Figure 28: The tensor network representation of the state \i]) v i) of Eq. [60 

We now look at the excitations in this model. There are both gauge and vertex excitations 
in this model. The vertex excitations are obtained by applying cr 2 and erf vertex v. The former 
commutes with the link operator and the plaquette operator but does not commute with the 
vertex operator A v creating a single vertex excitation and the latter commutes with the vertex 
and plaquette operator but does not commute with the four link operators adjoining a vertex 
v, thus creating four link excitations. Clearly these excitations are not anyonic like in the toric 
code case. The link excitations in fact depend on the valency of the vertices in the lattice. 
These excitations are shown in figure (I29jh 

The gauge excitations are obtained by applying erf and erf on the links of the lattice. The 
former commutes with the plaquette and link operators and creates charge excitations just as 
in the toric code case. They are deconfmed in the sense that applying a string of cr 2 operators 
on the direct lattice creates a pair of charges at the end points of the string irrespective of the 
size of the string. This feature is not present for the flux excitations as every time we apply 
a erf on a link we also create a link excitation on the link l. Thus if we separate a pair of 
fluxes by applying a string of erf operators on the dual lattice we end up creating a series of 
link excitations along the way. Thus the fluxes in this model are confined by the string tension 
provided by the link operator Ci. This is shown in figure (l30j) . 
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Figure 29: The vertex excitations in the matter sector of the H 2 /Z 2 model. 
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Figure 30: The gauge excitations in the H 2 /'L 2 model showing the deconhned charges and the 
confined fluxes. 


We now consider placing this model on a two dimensional lattice with a rough boundary. 
This changes the Hamiltonian at the boundary in such a way that they model is still exactly 
solvable. The bulk Hamiltonian remains the same and is given by Eq. [HHl The boundary 
Hamiltonian is given by 


H, 


rough boundary 


£ 


B r 


(61) 


pGdM 


with B p being the plaquette operators for the three sided plaquettes given by 


B p = 


cr 


a~ 


a 


(62) 


These operators are shown in the lattice in figure (T5TT) . It is clear that there cannot be any 
link operators on the rough boundary and hence the fluxes are now deconfined on the edge. If 
we try to move them into the bulk we have to overcome the energy barrier created by the link 
operator. Thus they are localized on the edges and are gapped. 

The existence of deconfined fluxes at the rough edge implies that the ground state degeneracy 
(GSD) of the system now increases. This is due to the loop operator 


L‘= ]J of (63) 

shown in figure (j52]h Thus the H 2 / r L 2 system when placed on a manifold with rough boundary 
develops a ground state degeneracy in contrast to the situation when it is placed on a closed 
manifold. 
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Figure 31: Deconfined fluxes localized at the rough boundary. 
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Figure 32: The loop operator acts on the non-degenerate ground state to give a new one. 

This feature is similar to the ones occurring in “confined Walker-Wang” models [31] as 
elaborated in [FT]. These models are exactly solvable 3D lattice models with a confined bulk 
but with deconfined anyons on the surface. Such surface states were also seen to be outside the 
group cohomology classification of SPT phases m in |43j 04]. The i/ 2 /Z 2 presented here is a 
2D model with this feature, the difference being that, only the fluxes are confined in the bulk 
and they get deconfined on rough boundaries. We believe these also to be outside the group 
cohomology classification scheme of SPT phases as the global symmetry does not act with an 
obstruction on these states [52]. Though we have seen this in this particular example it is easy 
to show that this is a property of all the systems we obtain from this formalism. In particular 
they are also true for non-Abelian groups. We will discuss this aspect of this model in another 
paper. 

5.2 ff 3 /Z 2 : 

The module H 3 is spanned by the basis elements {xoWTW}- The identity element of the 
gauge group Z 2 acts as identity on H 3 and is thus the three by three identity matrix. The other 
element 0_i flips x 0 an d Xi and leaves y 2 invariant. This makes the matrix representation of 
/i(0_ 1 ) as 

0 1 0 \ 

100. (64) 

0 0 1 / 
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Note that we can define other permutation actions of the gauge group Z 2 on H 3 where /i(0_i) 
leaves either xo or Xi invariant instead of y 2 ■ These actions are unitarily equivalent to the one 
defined above which leaves y 2 invariant. We will work with the one given in Eq. [64j 
The vertex operator in this model is given by 


Ay ~ 


1 + H v {(j)-i ) <g) a: 


cr. 


*2 


cr. 


*3 


cr. 


u 


(65) 


where a x is the Pauli matrix. The orthogonal operators become important in finding the 
excitations and hence we will also write them down. The operator orthogonal to Eq.fJHSD is 
given by 


At = 


1 /+(</* — l) 


at 


at 


a: 


*3 


at 


( 66 ) 


The link operator Ci is given by 


Q = 


1 + Z Vl + z 2 x 


1 + a; Z Vl + ujZ Vi 


1 + af 


3 

1 + Z Vl + z 2 x 

<8> 

3 

1 uj^ 1 Z V1 uj 

3 

1 + uZ vl + a J 2 Z 2 x 


1 - 


1 -of 


1 + af 


1 + Zy 2 + 


1 + Zy 2 + Z^ 


1 + UJ Zy 2 + U)Z V2 


1 + cu Z V2 + c oZ V2 


1 + LoZy 2 + UJ 2 Z% 2 


(67) 


where Z is a generator of Z 3 given by 


1 0 0 
Z= I 0 u 0 
0 0 CO 2 


( 68 ) 


and oj = e l tt. 

The plaquette operator is the same as the previous example. 

Once again we can write down exactly solvable Hamiltonians of the form 

« = (69) 

v p l 

with A v given by Eq.(j65]) and Ci given by Eq.(j67]) respectively. 

The ground states of this model satisfies the condition A v \gr) = B p \gr) = Ci\gr) = 
| gr) \/v,p,l. This is the same as in the previous example and is due to the fact that the 
Hamiltonian is made up of commuting projectors which leaves the spectrum of the model un¬ 
changed with respect to the previous example. Numerical computation of the trace of the 
transfer matrix gives us 5 as the ground state degeneracy for this model. We can write down 
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the 5 ground states of this model as 


IVb,o) = 

JI A v < 2 h |0i) < 8 >„ |xo) 

V 

(70) 

|Vb,2) = 

n Ay ®1 0 i) 1 x 2 ) 

V 

(71) 

1^,2 ,C*) = 

JI Ai ®iic{ | 0 i) ®i&ci | 0 -i) IX 2 ) 

V 

(72) 

1^,2 ,C*) = 

\A V 0 i) 0 -i) IX 2 ) 

V 

(73) 

\A ,2 ,c*,c*) = 

\Av ®l£C{,C% 0l) 0-l) X 2 ) 

V 

(74) 

where C{ and C 2 are non-contractible loops along the dual lattice around the two independent 
directions of the torus. It is easy to see that these states are linearly independent. Note that 
the hrst two states in this list have the gauge transformations acting in equivalence classes of 
the gauge group acting on the vector space carrying its representation, that is 7L 2 acting on 
H 3 . It can be seen from the gauge action that the equivalence classes in this case are given by 
{Xo 5 Xi} and {X 2 }- We can also see that within an equivalence class the state resembles that 
of the quantum double of the gauge group that is of the toric code in this case. 

We can write down another set of 5 ground states in the dual basis which are linear combi¬ 
nations of the above states. These are given by 

1 i) P ) = ik 

p 

Ci \ 4 >i + 4 >~t) ® v Ixo + xi + X 2 ) 

1 

(75) 

IV>A> = Y[b p 

p 

n^n +< ^- 1 ) ixo+ xi+ X 2 ) 

1 jec 1 

(76) 

10,^ 2 ) = l[B p 

p 

RCilla* ®i |0i + 0_i) ® v |xo + Xi + X 2 ) 

1 j&c 2 

(77) 

1 tCuCi) = n B r 

p 

Cl CTj ®1 01 + 0— 1 ) < 8 >„ Xo + Xi + X 2 ) 

1 jeCi,c 2 

(78) 

1 ifrvi) = IJA» 

I Ci ®i 10i) < 8 >„ Ixo + Xi + X 2 > 

(79) 


V l 


where C\ and C 2 are non-contractible curves along the direct lattice around the two independent 
directions of the torus. For example we can see that \ip p ) = + |0„,2) + \i/) v ,2,C*) + 

10,,, 2, (7|) + |0„, 2, C*, C 2 ). The other states obtained by applying the winding operator are 
similar combinations of the vertex states. 

Placing this system on a manifold with boundary gives new edge states, similar to the H 2 /'L 2 
case. The confined flux, in the {xo} equivalence class on the closed manifold, gets deconhned 
on a rough boundary. Apart from this deconhned flux we also have the deconhned hux in the 
{X 3 } equivalence class. However the ground state degeneracy decreases from 5 on the closed 
manifold to 3 on a manifold with a rough boundary. 
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5.3 tf 2 /Z 4 : 

We now consider another model with Z 4 gauge fields acting on the two dimensional vector space 
H 2 . The basis elements of Z 4 are denoted by {0o, 0i, 0 2 , 03} and those of H 2 are {xi,X-i}- 
The gauge action is given by the following two by two matrices, 


++) = h(02) 


1 

0 


0 

1 


and 

M0i) = M&O = 

This model has one extra deconfined flux in addition to the three deconfined charges. This can 
be seen due to the trivial action of the gauge group element 0 2 ° n H 2 . 

In terms of operators it is the statement that the operator Xf, where X is the shift operator 
given by 

/ 0 0 0 1 \ 

10 0 0 
_ 0 10 0 ’ 

\ 0 0 1 0 / 

commutes with the link operator Ci. This is seen to be true from the definition of the link 
operator. Consider 



1 Xvi 1 4*1 1 Xv2 ) 

(/^(^/) ’Xv\ \Xv2 )-^Z \Xv\ ? 0Z? XV 2 ) 

(80) 


= (Kb) -Xv i\xv 2 )\xv^ 4 * 1+2 , Xv 2 ) 

(81) 

whereas on the other hand we have 

\Xvi i&hXv-z) 

= (M+ 2 )^ \Xv 2 )\Xv^ 4>l+2, Xv 2 ) 

(82) 


(m(0z) 'Xvi \Xv2 ) \Xvi ) 4 *1+21 Xv2 ) 

(83) 


as 02 acts trivially on H 2 . Thus in this model we have deconfined an extra flux by a mere 
choice of representation space for the gauge group. 

5.4 Models with Deconfined Gauge Excitations 

The link operator Ci leads to confinement of the flux excitations in the previous two examples. 
Thus if we try to move a pair of fluxes apart we spend energy coming from the link excitations. 
The simplest way to avoid this is by making the link operator identity on all the links of the 
lattice, leaving only vertex and plaquette terms. The Hamiltonian for the H 2 /Z 2 case now 
becomes 

h = -Y j K-'£ b I < 84 ) 

V p 

where A* is given by Eq. fl42|) and B l p is given by Eq. flUj) respectively. This resembles the toric 
code Hamiltonian with the addition of the matter fields on the vertices. It is easy to see that 
the fluxes are no longer confined in this case as there is no cost in energy in moving a pair of 
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fluxes apart. The charge excitations are not confined as before. We have one vertex excitation 
as before on a vertex v which appears by acting with a z v on the ground state. 

We can carry out this construction for other gauge groups as well after choosing appropriate 
vector spaces for H n which carry the representation of these gauge groups. If we turn off the 
effect of the link operator Ci in all these cases we retain the deconfined flux excitations of the 
corresponding quantum double. 

5.5 Recovering the Quantum Double Hamiltonians of Kitaev 

If we choose the vector space on which the gauge group acts to be the trivial one dimensional 
vector space Hi, the gauge group acts trivially on such a space. This makes all the operators 
act non-trivially only on the gauge sector and effectively “switches off” the matter sector as the 
operators act trivially on this sector by just scalar multiplication. In particular the link operator 
Ci now becomes identity by construction. Thus we are left with only the vertex operator and 
the plaquette operators which act only on the gauge fields reproducing precisely the quantum 
double Hamiltonians of Kitaev. 


6 Outlook 

We have presented a systematic method to construct the transfer matrices of two dimensional 
lattice theories with gauge and matter fields inspired by the state sum construction of Kuper- 
berg [20]. The construction is a an extension of the one we started in ra by including matter 
fields. From the mathematical point of view the construction in H3 produced quantum doubles 
of various inputs of which gauge group algebras are special cases. In this paper we have more 
than a quantum double of these inputs due to the inclusion of the matter fields. If we use 
weak Hopf algebras and their modules as inputs in the construction presented in this paper we 
can obtain the Levin-Wen string-net models with confined excitations which was seen to be a 
feature in these models due to the presence of the link operator Ci. 

The exactly solvable models obtained for certain parameters in this paper were all sums of 
commuting projectors. They described fully interacting lattice theories with gauge and matter 
fields in two dimensions for arbitrary involutory Hopf algebras. We showed examples with long- 
ranged entangled ground states in this paper. We also expect to find models with short-ranged 
entangled states in this parameter space or by cleverly enlarging it to obtain the honeycomb 
model of Kitaev [46]. 

Moving from the ground state sector we find interesting excitations in these models. They 
contain both vertex and gauge excitations. While the vertex excitations are not anyons in 
the usual sense they nevertheless have interesting properties as was studied in [J5] where we 
showed that one can obtain non-Abelian fusion rules from a system with just Abelian degrees 
of freedom. While the charge excitations are deconfined as in the quantum double case, we 
have confined flux excitations in these models. The confinement arises due to the introduction 
of link operators which are symmetric under the local and global symmetries. Thus we have 
exactly solvable models with confined and deconfined particles with the latter protected by 
the energy gap created by these link operators. We will elaborate more on these models in a 
forthcoming publication. 
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A Input Data for the Construction 

The algebraic structures needed to construct the models comprise an involutory Hopf Algebra 
(-4, m, 77 , A, e, S) and a left n-dimensional ^.-module (. H , p) over a field C that is equipped with 
a co-structure t and a bilinear form G. We briefly describe each algebraic structure together 
with the properties relevant for the construction of the models. The notation we use is similar 
to the one used in m which was first introduced by Kuperberg in [20], 

A.l Hopf Algebras 

Let C be a field and A be an n-dimensional vector space over this field. Denote the basis by 
{0j}r= 1 an d its dual basis by { 0 *}” =1 , defined such that (fi l {(fij) = <5*. 

The vector field A is said to be an algebra if there are two linear maps: 

m : A £3 A — * A and p : C —> A, 

where m is an associative multiplication map and rj( 1) is the unit element. The multiplication 
map is defined by its action on the basis elements as follows, 

m{<fi a ® (fib) = m c ab (fi c , (85) 

where the sum over repeated indices is implied. The coefficients m c ab are called structure con¬ 
stants and can be thought of as tensors representing elements of A <8) A <8) A* so it is natural to 
associate them to the Kuperberg diagram on figure d33lh 

a 

m ab — 

b 

Figure 33: Kuperberg diagram for the multiplication map m. 

Also, we require the multiplication map to be associative such that the product of three 
basis elements (fi a , (fib and <fi c in A is: 


( ( fia,(fib)(fic (fia{,(fib(fic) j 

which can be expressed in terms of the structure constants as follows 

m k ab m l kc = m l ak m k bc: (86) 

Alternatively we can represent this associative property of the multiplication map using a 
Kuperberg diagram as shown in figure (l34jh 


a b a 



Figure 34: Associative property of the multiplication map. 


The unit element r] G A is such that for all x G A, xrj — rjx — x as depicted in figure (l35jh 
If both conditions shown in figures (l34jl and (1351) are fulfilled we say the triple {A, m, rj) forms 
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X 




Figure 35: Existence of the unit element for the multiplication map. 


an associative algebra with unit element r/. 

Similarly, we can take the dual vector space A* and define a multiplication map on it, A : 
A* ® A* — > A* by means of the action on two elements of the dual basis {0*}" =1 , 

A(0 4 ® <jP) = && = Aj>*, (87) 

where the coefficients A)( are the structure constants. There is an unit element that can be 
seen as a map e : C — > A* such that e(l) is the unity for the dual multiplicative map A. So 
the triple (.A*, A,e) defines an algebra structure in A *; Equivalently we can regard the maps 
A and e as being a co-multiplication and co-unity in A respectively. In this sense, the map 
A : A —* A <S> A is defined as 

A(</>*) = A® 0 fc ), (88) 

and the unit map e : A —> C is defined as 

e(&) = eh (89) 

The triple (A, A, e) forms a co-associative co-algebra with co-unit provided the following two 
relations are fulfilled 

k i > k > A 

a b a 

(a) 

Figure 36: Associativity of the co-multiplication map and existence of the co-unit of A. 

The quintet (A, m , 77, A, e) is said to form a bi-algebra whenever some special compatibility 
conditions are satisfied. This is, when the co-multiplication and the co-unit map are homomor- 
phisrns of the algebra, namely 



A {(j} a ^b) = A(0 a )A(0 fe ), (90) 

e{4>a4>b) = e(0a)e(0&). (91) 

The Kuperberg diagrams for these conditions are shown in figure fl37jl . 

Consider now the endomorphism S : A —> A called the antipode of the algebra. If such a 
map satisfies the condition shown in figure (1351) then we say the sextet (A, m, r], A, e, S) forms 
a Hopf Algebra. Additionally the sextet can be equipped with a ^-structure which ensures the 
definition of Hilbert spaces over complex numbers and unitarity of the system mm- Consider 
the conjugate linear involution * : A —> A satisfying 

(x*)* = x, 

(xy)* = y*x*, 

V* = V 
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(92) 

(93) 

(94) 



(a) (b) 



A^pm 

A : -Am 


( c ) 


Figure 37: Bi-algebra compatibility conditions. 




Figure 38: The antipode condition. 


for all x,y G A and the unit element rj G A. Then the triplet (A, m, rj) is called a ^-algebra. 
Moreover, the quintet (A, m, A, r/, e) is called a *-bialgebra if the involution map * is consistent 
with the co-product structure, i.e. 


A(x*) = A(VP, (95) 

e(x*) = e(x), (96) 

where the bar in the r.li.s. of the last expression stands for the complex conjugate. Finally, 
the sextet (A, m, A, 77 , e, S) is called a Hopf *-algebra if the conjugate linear involution * is 
compatible with the antipode S in the following sense: 

S(S(x)*)* = x. (97) 


A.2 Left A-Module 

Consider the triple (A, m, rj) as being an algebra over the field C, let H n also be a vector space 
over the same field C, and let y : A <S> H n —>■ H n be a morphism of vector spaces such that 
it takes an element g G A and acts with it on an element y G H n , we denote this action as 
/i( g ) > x- The pair (H n , g) is said to be a left A—module pfE, Mj if the following properties are 
satisfied for all y, y' G H n \ c G IK, g, h G A and r) the unit element of the algebra. 


MsO e H n , (98) 

H(hg) > y = g(h) > (g(g) > y), (99) 

V>(v)>X = X, (! 00 ) 

M > (ex) = C (y(g) > x) ( 101 ) 

M > (x + x') = M > x + Kg) > x' (102) 


Let us associate the diagram in figure (j5U|) to the new structure constant of the left CA—Module. 
Note that we make a distinction between the kind of arrows in figure (1391) . This is due to the 
fact that the module map /i combines elements of two different spaces. Consequently we use a 
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Figure 39: Diagram associated to the module structure. 


green arrow to represent an element of the vector space H n while a black arrow represents an 
element of the vector space A. 

Using this diagrammatic representation we can re-write properties in Eqs. (1981) to (11021) as 
depicted in figure (l40l) . 





(a) Ea. IM1 



u+x'y 


-JJ- 


(b) En.rmnl 




(c) Ea. IlHlI 


(d) Ea. 11021 


Figure 40: Diagrammatic notation of properties IfFJl to 11021 


We equip the pair (H n , /i) with a co-associative co-multiplication structure given by the 
map t : H n —)■ H n ® H n . If {Xa}a=i are the basis elements of H n the co-multiplication map is 
dehned by the action on these elements as follows: 

t{Xa) = X 7 ), 

where are the structure constants and are represented by the Kuperberg diagram in figure 
(141( a)). Moreover, the co-structure is co-associative meaning that the condition depicted in 
figure (TUT b)) is satisfied. 



(a) (b) 

Figure 41: (a) The co-structure tensor, (b) Associativity of the co-product 

There is a compatibility relation between the left M-modulc and the co-product of the Hopf 
algebra A shown in figure (1421) . It implies that the co-structure map t preserves the action of 
the A- module, i.e., 

>x) = M A (°)) >t(x) (103) 

Finally, we define a bilinear map G : H n <g) H n —y IK. It can be represented by a tensor 
with two incoming green arrows as shown in figure (1431) . This tensor will help in contracting 
the arrows for the matter degrees of freedom when constructing the partition function of the 
model. 

The above algebraic structures are the elementary building blocks of the partition function 
and the transfer matrix from which we will obtain specific models written in terms of local 
operators. 
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Figure 42: The compatibility relation between the left A-module and the co-product of the 
Hopf algebra. 

G a p = a -- P 


Figure 43: The bilinear map G and its tensor representation. 


B Splitting of the (1 + 1)D Transfer Matrix 

B.l Associating Tensors to the Lattice 

In this case we consider a 2-dimensional lattice C composed of vertices, links and faces as shown 
in figure (144]) 






• - • 

Figure 44: A (1 + 1 )D square lattice cell and its components, vertices, links and face are shown. 

As in § (1.1.1) to each of the lattice components we associate a set of tensors. The M aia2a3a4 
tensor is associated to the faces of the lattice. Since each vertex C is now four-valent a 4-tensor 
T q i Q 2Q3Q 4 j s associated to it. To the links of the lattice we associate a tensor, as each link 
is connected to two vertices and two faces. This is schematized in figure (145|h 

The M aia2a3a4 tensor involves gauge degrees of freedom only and consequently has Latin 
indices only. The T ai0l20l3aA involves only matter degrees of freedom therefore it has only Greek 
indices. Consequently, the mixed tensor involves contractions between matter and gauge 
degrees of freedom. The contraction of these tensors is orientation dependent. The antipode 
of the Hopf algebra S : A —* A takes care of the orientation in the gauge sector whereas the 
bilinear form G : H ® H — > C performs the task for the matter sector. Thus, the contraction 
rules are the same as in the (2 + 1 )D case. For the gauge degrees of freedom an antipode is 
placed whenever the orientation of the plaquette does not match the one of the link. Similarly 
for the matter degrees of freedom, the green arrow coming out from the mixed tensor is 
joined to the T a i"2«3«4 v j a the bilinear form Gp ai . 

B.2 Splitting the Transfer Matrix 

The splitting procedure shown in 13.31 can, as well, be used to write the transfer matrix of the 
(1 + 1 )D case as a product of local operators. The result of contracting all the tensors on the 


42 













= “1- *-M~* -“3 


(a) 


a 4 

OL 1 - • - a 3 = rj-iOLlQ'2 0t3&4 = a>i~* - 'J 1 - ^ 0-2 

OC2 


(b) 


= L abf3 a = 


a ■* - L -'■# <> 


(c) 


Figure 45: The tensors associated to each component of the lattice. 


lattice is defined to be the partition function of the theory, given by: 


Z[£,i A., Hi a>i 11 -^aia2<i3a4 ip) n rjiai 


OL20L3CX-4 ( 


1^(0 IMG 


a/3* 


(104) 


tk 


z= 




i A A i A 


Figure 46: The graphical representation of the partition function Z(£,A,H,z,£) is shown. 
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The partition function is related to the one step evolution operator U as follows: 


Z = tr(U N ), 


(105) 


therefore U is represented as a single slice on the graphical representation of figure f!46[) . 



Figure 47: The one step evolution operator U(C,A ,) is one slice of the partition function. 

From the correspondence shown in figure (l45|) the one step evolution operator can be written 
as a tensor network built from the tensors representing the structure constants of the algebra 
A , the module 77, the antipode S and the bilinear form G. This results in a rather intricate 
representation that will be split in operators acting on smaller pieces of the total Hilbert space, 
as we will we show later on. 
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Figure 48: The tensor network representation of the one step evolution operator is shown. 

The process of writing the one step evolution operator U as a product of local operators is 
based on the algebraic properties of the tensors that represent the elements of the lattice. We 
begin by splitting the diagram in figure (147)) into two parts which we call A and B as shown 
in figure (1491) . Now we proceed to the splitting of the A part. It is sufficient to pick one of 
the plaquettes and the two adjacent timelike links and write the tensor network representation. 
This is depicted in figure (1501) . 

The associative property of the multiplication map m allows us to literally split the M 
tensor representing the plaquette of the lattice as in figure (l5lT a)). At this point, we are able 
to define the vertex operator as the tensor network shown in figure (151( b)) such that it acts on 
two links denoted A, i 3 and the vertex i 2 in between. 

Therefore, the U operator now can be represented as a product of vertex operators A v for 
all vertices v in the lattice, as depicted in figure (152|) 
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Figure 49: Step 1 of the splitting procedure 



Figure 50: Tensor network representation of a plaquette and the two adjacent timelike links. 
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Figure 51: (a) Splitting of the M tensor, (b) The vertex operator A v 


Similarly, we now proceed to the splitting of the upper section in figure (149]1 that we called 
B. It is enough to consider two vertices and the spacelike link in between, consequently the TN 
representation is shown in figure (T53h . 

The associative nature of the co-multiplicative map t allows us to split the T tensors on 
each vertex as shown in the l.li.s of figure (153f b)) by writing the mixed operator in terms 
of its elementary components we get the link operator C) as the TN highlighted in the r.h.s. 
of figure (I55f b)). Thus, the B part of the U operator is now written as a product of C) for 
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Figure 52: The U operator after the splitting of the A part of the tensor network in figure (142)1 . 
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Figure 53: In (a) we show the tensor network representation of a pair of vertices and the 
spacelike link in between. In (b) the link operator is explicitely shown. 


each link of the lattice together with the connector operators which we call Li and V v that act 
on single spacelike links and vertices of the lattice, respectively. Finally the one step evolution 
operator has been decomposed as a product of the local operators depicted in figure (15T) . 
Therefore, the complete transfer matrix of the theory can be written as: 

U (zt, Zg, Zt, my, Gs , Gt) = Ci(Gs) n v n >■ n A v (z* t ) (106) 

l l V V 

where the tensor network representation of each local operator is shown in figure (l55j) . The 
operators V v and Li can still be written as shown in figure (|56j) . in which we have used the 
operators defined in figure (1211)1 and the antipode identity for involutory Hopf algebra. Finally 
we can write the transfer matrix as 

U An, Zgi Zrpi tTL V i Gs') 11 Ci(Gs) Qv{j^v) U(zT)r,(4) IW Z W ■ 

l V l V 

The operator Di can shown to be proportional to identity, more precisely Di = |G|I, while the 
operator V v (Gt ) can be made into identity by choose (Gt)^ 3 — 5(a,(3). 
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Figure 54: Final form of the U operator after the systematic split of the TN performed over 
the entire lattice. 
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Figure 55: The tensorial representations of each local operator that make up the one step 
evolution operator U. 


C Examples of Exactly solvable 1-D Quantum Models 

In these examples there are no plaquette operators as there are no plaquettes in one dimension. 
The models we consider will be made of only the vertex operator A v , which are again the 
gauge transformations, and the link operators Ci which describe the gauge and the matter 
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(b) operator which acts on a single link. 

Figure 56: Operators which act on the links and vertex of the lattice. 


interactions. Thus our Hamiltonians will be of the form 

h = -J2 a v-J 2 Ci - ( lor ) 

V l 

We will just look at two examples, // 2 /Z 2 and F/ 2 /Z 4 . We will see that in the first case there 
are deconhned vertex and flux excitations with no ground state degeneracy. In the second case 
there continues to be deconfined vertex excitations but the ground state in this case has a lot 
of degeneracy which we attribute to the “condensation” of certain fluxes of the corresponding 
model in two dimensions. 

C.l H 2 /Z 2 : 

The action of 7L 2 on H 2 is the same as in the two dimensional example. 

The vertex operator is given by 


A 


V 


1 + 


(108) 


where l\ and I 2 are the two links adjacent to the vertex v. 
The link operator is given by 


Ci = 


1 + 


z 

V2 


2 


(109) 


where V\ and V 2 are the two vertices flanking the link l. 

It is easy to see that these operators commute with each other and are projectors. Thus 
as before the Hamiltonian is exactly solvable and its spectrum is easily obtained, the ground 
state conditions are similar to the corresponding example in the two dimensional case which is 
A v \gr) = Ci\gr) = | gr), Vw,Z. Using the arguments in the two dimensional case the number of 
ground states is again one for this model. They are constructed in manner similar to the way 
it was done for the two dimensional case. They are given by 


IVb) 

= JJa, ®i \<t>i) Ixi) 

V 

(110) 

\^l) 

= Y[Ci ®i 101 + 0_i) <g>„ \xi + X-i)- 

(111) 
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As before we expect these two states to be the same. 

The excitations are obtained as before by applying a x to either the vertices or the links 
to create link excitations. The former gives us deconhned link excitations and the latter gives 
us isolated link excitations. By applying a z to either vertices of links we obtain the vertex 
excitations. The former gives us isolated vertex excitation whereas the latter gives us deconhned 
charge excitations. The isolated vertex and link excitations are shown in figure (17T|) and the 
deconhned gauge excitations are shown in hgure (l58|i respectively. 

a v a l 

- * -•-•-•— X —•-•- ••• 

Figure 57: The isolated vertex and link excitations in the model. 


Figure 58: The deconhned vertex and link excitations in the H 2 QL 2 model. 


C.2 ^ 2 /^ 4 * 

The gauge action of Z 4 on H 2 is taken to be the same as the one in the two dimensional 
case. The novelty in one dimension is an interesting phenomena by which the ground state 
degeneracy increases. This is because a single operator X? on a link l commutes with the 
Hamiltionian. This was shown in the corresponding example in the two dimensional case, the 
difference here is that there is no plaquette operator in one dimension. Note that this is not 
a gauge transformation as then we would have to act on the adjacent vertex and link of l as 
well. Thus we have effectively “condensed” the two dimensional fluxes obtained by applying 
the operator Xf along a path in the dual lattice. The deconhned charge excitations still exist 
thus giving us a model with deconhned excitations and ground state degeneracy more than one 
in one dimensions. 
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